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Abstract The dipole-dipole interaction influences nanoscopic matter by fixing the patterns 



of permanent, displacive, and induced dipole moments, subject to constraints of 
molecular size and other short range interactions. Prediction of these arrange- 
ments is a challenging problem. The eigenvector of maximum eigenvalue of the 
dipole-dipole interaction matrix can provide insights and sometimes a complete 
solution. As an example, the octahedral tilt instabilities of perovskite-type crys- 
tals is shown to optimize dipolar interactions. Therefore this instability can be 
designated as antiferroelectric, because dipole-dipole interactions are a domi- 
nant driving force. 
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1. Introduction 

In condensed matter one finds three kinds of electrical dipoles: (1) per- 
manent (as in a molecular cluster of polar molecules, such as CH3CI), (2) 
diplacive (as in BaTiOa), and (3) induced (a cluster of benzene molecules in 
an external field). The unifying feature is that the dipole moments feel each 
other. Their magnitudes and orientations depend on where the other dipoles 
are, in a self-consistent fashion. This paper argues that one should study the 
linear algebra, and particularly the eigenvector of maximum eigenvalue, of the 
dipole-dipole matrix V which relates the induced electric field -F^md at the i'th 
site Ri to the dipoles fij at sites Rj : 
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where Rij = R4 — Rj. The insights obtained from the eigenvectors and eigen- 
values are illustrated below in some simple examples, and in an application to 
real materials with spontaneous electrical polarity, the octahedral tilt instability 
of Re03-type perovskite crystals. 

Of course, it must be acknowledged that not all properties of polar matter 
have much to do with the dipole moments of their molecules. For example, the 
water molecule has a dipole in vapor phase, and electrical polarity is a dom- 
inant effect in the interactions of water molecules with each other and with 
dissolved species. However, for near-neighbor pairs, it is not a good approxi- 
mation to replace the polar charge distribution of water by a single dipole. A 
better approximation would use two dipoles, located along the two bonds, each 
pointing away from the oxygen toward one of the two hydrogens. The sum of 
these two dipoles is the total measured dipole moment of the water molecule. 
Also, it would be not a very good approximation to assume that the two dipole 
moments of the water molecule are unchanged in their various chemical envi- 
ronments. Induced polarity should be treated on top of permanent or displacive 
polarity, and this is done in the case of the octahedral tilt. A complete theory 
would not make any type of dipolar approximation, but would deal with the 
actual varying charge distribution of the molecules in detail. Modern density 
functional theory does this well, but slowly, which motivates a search for a less 
complete but still sensible theory. 

2. One Polarizable Molecule 

The energy of a single molecule is n 2 /2a — fl ■ F cxt where a is the polariz- 
ability and F cxt is an externally applied field. The actual dipole moment jl is 
fixed by the condition that the energy is minimum, which gives fl = aF ext . 

3. More than One Polarizable Molecule 

Now consider N molecules, approximated as points, at chosen fixed posi- 
tions Rj,. It costs energy fJ.f/2a to create a dipole on the ith site, but one 
gets back energy —jXi ■ Fi from the total field F{ = F i in ^ + F, iCxt at the site, 
where F{ )Cxt may be spatially inhomogeneous, varying from site to site. The 
total energy, in a standard vector space notation, is 
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When there are N sites, each with a dipole, then \fj, > is a 3N column vector 
containing all the components fii a for a = x,y,z and i = 1, . . . , N. Notice 
that the interaction energy — pi ■ F^md is reduced by a factor of 2 in the total 
energy, to avoid double counting. With no external field, the system is stable 
if 1/a — T is a positive operator. Equivalently, all eigenvalues of T should be 
less than 1/a. 

Just as for the single molecule, the induced dipoles pi should minimize the 
total energy, dU /dfii a = 0. The solution is 

1 ^ 



| A* >= ( - " "J l^ext > (4) 

This is a generalized version of the Clausius-Mossotti law [1-3]. An alternate 
notation, |/i >= XXV a ~~ 7) < Tl-^ext > |7 > uses an expansion in 
eigenvectors, where I7 > is the eigenvector of T with eigenvalue 7. If the 
overlap < 7max|-Fcxt > is not zero, the eigenvector of largest eigenvalue 7 max 
may dominate the response. 

If 7max increases to become equal to 1/a, then the system becomes un- 
stable. It will develop a spontaneous polarization, with a pattern |/i > pro- 
portional to the corresponding eigenvector |7 max >■ Higher order terms not 
considered here are needed to determine the magnitude of the spontaneous 
moment. 

4. Two Dipoles 

For pedagogical purposes, let us examine the 6x6 matrix V which describes 
the interaction between two dipoles separated by distance \R\ = a; to be spe- 
cific, let R point along the z axis. Even without writing the matrix, we can 
find all six eigenvectors by simple physical reasoning. The eigenstates of V are 
those patterns of dipoles such that the induced field at each dipole is parallel 
and proportional to the dipole itself. Consider the patterns in Fig. 1 .a-d. The 
field F\ lies along the z-axis if p2 points along z, and lies in the xy plane 
antiparallel to p2 if P2 lies in the xy plane. Therefore vectors lying along z 
separate from those lying in the xy plane, and by cylindrical symmetry, eigen- 
vectors in the xy plane are doubly degenerate. The degenerate pair can be 
chosen so that one lies along x and one along y. To save space, we project the 
problem onto the 2d xz plane. The resulting 4x4 matrix is 
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7^ 
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The eigenvalues are (2,1,-1,-2) in units of (1/a 3 ), and the corresponding eigen- 
vectors, shown in Fig. l.a-d, are 



(6) 



These eigenvectors all have different symmetries. Under inversion and mirror 

reflection in the central xy plane, the vectors have (Aj, A m ) = ( ), (H — ), 

( — h), and (++) symmetry, respectively. 
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Figure 1. Dipole patterns of two and three interacting dipoles. In a-d the eigenvectors of V are 
shown for two dipoles. In e,f,h, the eigenvector of largest eigenvalue is shown for three dipoles. 
In e, the dipoles all have equal length. In g and i the optimal patterns (not eigenvectors!) for 
fixed length dipoles are shown, nearly agreeing with the eigenvectors of f and h. The vectors in 
h and i differ in angle by about 3° . 

The isotropic Heisenberg interaction U = J pi • \ti for two spins has exactly 
the same eigenvectors but very different eigenvalues. For J < (ferromagnetic 
exchange), states a and c are low energy states, equal in energy, and states b 
and d are high energy states, also equal in energy. For J < (antiferromag- 
netic exchange), the reverse is true. For magnetic systems (the O2 molecule 
is a nice example), the exchange energy is much larger than the dipole-dipole 
interaction. The ground state of O2 has two outer electrons, polarized ferro- 
magnetically (S = 1), and only a tiny spin-orbit-induced preference for the 
moment to point perpendicular to the axis. By contrast, electrical dipoles have 
a large dipole-dipole interaction which is intrinsically anisotropic and strongly 
prefers state a. 
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5. Three Dipoles 

In the two-dipole problem of the previous section, the eigenstates all had 
the property that each molecule's moment /2j (where /ij a =< ia[y >) had 
the same magnitude |/2j|. This brings a nice benefit: these patterns tell us also 
about the energetics of the problem of fixed size permanent dipoles. 

Fig. l.e-i shows geometries and corresponding dipole patterns for three 
molecules. The physical situations of interest are (1) polarizable molecules 
with no permanent moment, (2) permanent dipoles which can orient in space 
but not change their magnitude, and (3) permanent dipoles which are also po- 
larizable and change their magnitude. We omit category (3) for now. Category 
(1) is mathematically simplest because linear algebra determines what happens 
in an external field, and whether spontaneous polarity appears. Fig. l.e,f,h 
show the eigenvectors of largest eigenvalue of V for (e) the equilateral triangle, 
(f) the linear arrangement, and (h) an obtuse isosceles triangle. For the equi- 
lateral case, the eigenvector is determined by symmetry, but for the other two 
cases, there are two independent vectors odd under the perpendicular mirror, 
so a 2 x 2 matrix is needed, and the central atom (or molecule) has a larger 
|/2i n d| than the end atoms. 

Because \ j2i\ is no longer the same for each molecule in the maximal eigen- 
vector, it is not the solution of the problem of a fixed moment. Instead, for 
fixed moments we want to know the pattern \fi > which minimizes 

U(fli, . . .,fl N ) = ~\< ^l r lM > - < n\F ext >, (7) 

subject to N nonlinear constraints \ = fiQ. This is a much harder mathe- 
matical problem. Fortunately, the solution is guaranteed [4] to have a smooth 
relation to the maximal eigenvector of l~. The solutions to the fixed moment 
problem in zero field are shown in Fig. 1 .g,i, and are close to the unconstrained 
solutions in Fig. l.f,h. 

6. Infinite Stack of Dipoles 

Consider the arrangement \fi > of Fig. 2, with equally spaced dipoles /2of 
equal magnitude, all pointing in the z direction which is also the direction of 
the stack of molecules. It is clear that this pattern is the maximal eigenvector 
of the corresponding l~. The eigenvalue is 4£(3)/a 3 where a is the spacing 
of dipoles and £(3) = 1.202057 ... is the Riemann zeta function. Since the 
moments are equal, this pattern is also the solution to the problem of fixed 
dipoles on a 1-d line. It is interesting to compute the electric field of this stack, 
which is the negative gradient of the electrostatic potential <J>(p, z) (where p = 
\Jx 2 + y 2 is the radial distance). By methods described elsewhere [4] the 
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Figure 2. Faraday flux lines indicating the field strength and direction (indicated by arrows) 
for a single dipole pointing up (left panel) and for a periodic stack of dipoles pointing up (right 
panel). In the left panel the field falls as 1/r 3 in all directions. In the right, the field decays 
exponentially in the directions transverse to the stack. 

potential has the asymptotic form 



This result explains why fixed dipoles particularly like to form one-dimensional 
stacks. The dipole-dipole attraction within the stack is very big, and the inter- 
action with other dipoles falls very rapidly with distance. The asymptotic form 
Eq.8 is already accurate to 2% at p = a. 

7. Cubic lattice of Dipoles 

An eigenvector of this problem was found independently by Clausius [1] 
and Mossotti [2], and is illustrated in Fig. 3. If all dipoles point along a cube 
axis, with equal magnitude, it is clear that the induced field does the same 
thing, so this is an eigenvector. The Clausius-Mossotti arguments show that 
the eigenvalue is 7cm = 47m / 3 where n, the density of molecules, is 1/a 3 
for the simple cubic lattice. A good pedagogical discussion is given by Aspnes 
[3]. One should also ask whether this is the maximal eigenvalue. From Fig. 
2, it is clear that for the simple cubic lattice illustrated in Fig. 3, the answer is 
"no." One can further lower the energy by reversing the direction of alternate 
stacks of dipoles, producing a slightly larger eigenvalue in an antiferroelec- 
tric pattern. Luttinger and Tisza [5] proved that this pattern had the extremal 
eigenvalue. For diamond structure, a more complicated extremal antiferro- 




(8) 



Bled - Dipoles 



1 



Figure 3. Infinite periodic lattice of cubic symmetry, showing the ferroelectric eigenstate 
studied by Clausius and Mossotti. Also the 2-d analog is an eigenstate. These are not necessarily 
the maximal eigenstates. 



electric eigenstate was found by White et al. [6]. These authors also showed 
that for fee and bec structures, the Clausius-Mossotti ferroelectric eigenstate is 
maximal. 



8. Octahedral Tilt 
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Figure 4. Energy versus tilt angle for (7r,7r,7r) periodic tilts around the (111) axis of BX3 
(ReC>3 structure type) perovskites. The solid line is the total electrostatic energy (units Z 2 e 2 / a, 
where Z is the anion charge) computed for rigid ions comprising rigid octahedra, with a corre- 
sponding hard-core contraction perpendicular to (111) planes. The dashed curves , from top to 
bottom, are electrostatic energy for one displaced anion, dipole-dipole attraction energy of in- 
teracting displacements, and volume contraction energy. The sum of these is the dashed curve, 
coinciding with the exact energy up to fourth order terms in the Taylor series. 



8 



In this last section a brief discussion is given of recent work to be explained 
more fully elsewhere [7-9]. We have found that dipole-dipole interactions 
play an unexpectedly big role in the "octahedral tilt" instabilities common in 
perovskite materials [10-13]. The motif is evident in Fig. l.e, showing the ten- 
dency of a triangular plaquette to polarize circumferentially. In the perovskite 
tilt, when the axis is (111), the displaced anions create the same triangular pat- 
tern of displacive dipoles. The moment points from the new anion position 
to the undistorted cubic anion position. These moments interact because they 
create an induced field at each undisplaced site. After some thought one can 
conclude that this is an eigenstate of l~, and after some calculation one can 
prove that it is the maximal eigenstate. There is no better way to place dipoles 
on the X sites of ABX3 perovskite than in the pattern corresponding to this tilt. 
The eigenvalue 7 max = 14.461/a 3 , is 15% bigger than 7cm- The consequence 
is that in a rigid ion picture, the BX3 simplification of the ABX3 lattice is only 
marginally stable under electrostatic and hard sphere forces. When the anion 
polarizability is added to the theory [7], the structure is almost immediately 
unstable even for very small polarizability. The maximum eigenvalue is acci- 
dentally 5-fold degenerate. Three of the partners permit (tt, tt, tt) alternating 
tilts in any linear combination around any cube axis. These generate the dis- 
tortions labeled a~6~c~ in Glazer notation [10]. The Glazer distortions with + 
signs require tilts with (71", tt, 0) periodicity, which happen also to be eigenvec- 
tors of T, with slightly smaller than maximal eigenvalue. It seems appropriate, 
given the role that dipolar interactions play, to regard all these distorted ground 
states as multiple-sublattice antiferroelectrics. The other two partners are re- 
lated to the cooperative Jahn-Teller ground state of LaMn03, which thus also 
has an unexpected contribution from polar interactions in its distorted ground 
state. 

9. Summary 

Maximal eigenstates of the dipole-dipole interaction tensor V are useful and 
relevant to several problems: response of polarizable clusters to external fields, 
causes and patterns of spontaneous distortion of polar assemblies, and the self- 
organization of polar molecules in various media. It is surprising how robust 
and useful this concept seems to be. Given that short range repulsive forces 
are large, and that near-neighbor electrostatics may involve higher multipoles, 
there is no obvious reason why a pattern like the perovskite octahedral tilt 
should have to conform to the maximal eigenvector. The fact that it does should 
encourage us to use this tool in other problems. 
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